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Abstract. We study the relativistic dynamics of a pressure-less and irrotational fluid of 
dark matter (CDM) with a cosmological constant (A), up to second order in cosmological 
perturbation theory. In our analysis we also account for vector and tensor perturbations and 
include primordial non-Gaussianity. We consider three gauges: the synchronous-comoving 
gauge, the Poisson gauge and the total matter gauge, where the first is the unique relativistic 
Lagrangian frame of reference, and the latters are convenient gauge choices for Eulerian 
frames. Our starting point is the metric and fluid variables in the Poisson gauge up to second 
order. We then perform the gauge transformations to the synchronous-comoving gauge and 
subsequently to the total matter gauge. Our expressions for the metrics, densities, velocities, 
and the gauge generators are novel and coincide with known results in the limit of a vanishing 
cosmological constant. 
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1 Introduction 

To put it in simple words, cosmological structure formation is nothing but fluid dynamics. 
There are two frameworks to study fluid dynamics, namely the Lagrangian and the Eulerian 
frame. The Eulerian frame of reference is associated with an observer’s coordinate system, 
where the observer is studying how the distant streams of matter are clustering. The La¬ 
grangian frame of reference, by contrast, is the physical frame attached to the clustering 
matter elements, i.e., each observer is literally following an individual matter element. In the 
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simple Newtonian world, there is only one Eulerian frame, because the infinite Eulerian co¬ 
ordinate systems are physically equivalent, thanks to the Galilean invariance of the Eulerian 
fluid equations. In this Newtonian set up, the Eulerian and Lagrangian frame are uniquely 
defined by the spatial coordinate transformation embedded in an Euclidean geometry 

x* = (? + &, ( 1 . 1 ) 

where the bar denotes Euclidean coordinates, x l {f,q) is the Eulerian coordinate, q 1 denotes 
the Lagrangian coordinate (i.e., initial position) of a given fluid element, and is the 

Lagrangian displacement field. 

In General Relativity (GR), however, this situation is more involved, mainly because of 
the covariant formulation of the relativistic equations of motion. To be more specific, the 
fully relativistic Lagrangian frame can be uniquely identified with the synchronous-comoving 
gauge, because the resulting space-time coordinates are indeed the ones attached to the 
matter [1], But there are infinite possibilities for selecting an Eulerian frame in GR. Even 
defining an Eulerian frame in GR is not so straightforward. Naturally, perhaps because our 
thinking is Newtonian, one could define a generic Eulerian frame in GR with the following 
correspondence to the Newtonian theory: 

In General Relativity, an Eulerian frame can be identified with a certain gauge if 

the resulting equations of motion take the classical Newtonian form. 

As a direct consequence of this definition, it follows that the gauge generator of the gauge 
transformation 

x u = q u + ^ (L2) 

is the 4D relativistic counterpart of the 3D Newtonian displacement field (cf. with eq. (1.1)). 
Here, are the Eulerian coordinates of a given Eulerian gauge, and q^ = {r,q l ) are 

the coordinates of the synchronous-comoving gauge. With the above quote, we fix the spatial 
coordinates of a generic Eulerian frame. To get a specific Eulerian frame in GR, one then 
needs to fix the time coordinate, which is in principle arbitrary. One particular convenient 
choice of the time coordinate is certainly to use the proper time of the particles (which is 
precisely the case in the total matter gauge, see section 4). 

All gauges in GR are physically equivalent, but the physical measurement of the density 
and the velocity of the matter streams depends on the particular Eulerian coordinate system 
chosen by a given observer. For example, convenient Eulerian gauges are the Poisson gauge 
and the total matter gauge, the very Eulerian gauges we shall investigate. The Poisson gauge is 
particularly common for studying the gravitational lensing [4], whereas the total matter gauge 
is important for the Eulerian biasing, and when relating GR calculations to investigations in 
a flat space-time, such as Newtonian IV-body simulations [5-7]. The synchronous-comoving 
gauge, apart from being the unique Lagrangian frame of relativistic fluid flow, is the gauge 
for studying the Lagrangian biasing, especially when considering the effect of the primordial 
non-Gaussianity [8-13]. 

In this paper we consider the cosmological fluid in the above mentioned gauges in the 
framework of cosmological perturbation theory. (For different solution techniques, see e.g., 
refs. [14-21].) In detail, we study the relativistic dynamics of a pressure-less and irrotational 

1 The locally inertial (or Minkowskian) Lagrangian frame, by contrast, is represented by the locally free- 
falling coordinate system along the worldline of the observer conroving with the particle, i.e., the Fermi frame, 
see e.g. [2, 3]. 
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fluid of dark matter (CDM) with a cosmological constant (A), up to second order in cosmo¬ 
logical perturbation theory. We begin with reporting the results of refs. [11, 22] in the Poisson 
gauge, simplify their expressions and provide explicit solutions for all time coefficients, and 
then we obtain the results in the synchronous-comoving gauge and in the total matter gauge 
by means of gauge transformations. We generalise to a ACDM Universe with primordial 
non-Gaussianity the results for the synchronous-comoving and the Poisson gauge of ref. [23], 
which were restricted to an Einstein-de Sitter (EdS) model, and extend the analysis to the 
total matter gauge. In this paper, almost all of the results are new, only few results were ob¬ 
tained before in the literature. For example, the matter density in the synchronous-comoving 
gauge was obtained in refs. [6, 11, 24-26], whereas the matter density in the total matter 
gauge was obtained in refs. [6, 24-26]. In the literature, there are no known ACDM results 
for the second-order metric in the synchronous-comoving and for the total matter gauge, ex¬ 
cept for ref. [6] for scalar and vector modes, however obtained with the gradient expansion 
technique. Here we find the full metric results, the densities and the velocities for these three 
gauges. In our analysis, we pay particularly attention to distinguish between Newtonian-like 
and GR-like contributions. Note that the fluid approximation and thus our approach is only 
valid on sufficient large scales, i.e., before matter streams start to intersect (“shell-crossing”). 
Furthermore, our analysis is restricted to the fastest growing mode solutions, i.e., we discard 
decaying modes. 

This paper is organised as follows. In the following section we introduce our conventions 
of the metric, density and velocity perturbations. In section 3 we give the general formulas 
for gauge transformations up to second order. Then in 4 we give the definitions of the used 
gauges. The derivation of our results in the Poisson gauge are given in section 5. This analysis 
serves as our basis to perform the gauge transformation to the synchronous-comoving gauge 
(section 6). In section 7 we then transform the Lagrangian results to the total matter gauge. 
The former two sections are intentionally kept quite pedagogical. For the readers interested 
only in the results we refer to the compact section 8. Finally, we conclude in 9. 

Notation: we use /i, v... for space-time indices and i,j ,... for spatial indices. Summation 
over repeated indices is assumed. If not otherwise stated, spatial indices are raised and lowered 
with the Kronecker delta. A comma denotes spatial partial derivative and a dot denotes partial 
derivative with respect to the conformal time. A bar denotes from here on FLRW background 
quantities. We set c = 1. Quantities in the synchronous-comoving gauge, in the Poisson 
gauge, and in the total matter gauge are respectively denoted with the subscripts S, P, and T. 
Consistently, quantities which arise because of a gauge transformation are abbreviated with 
the respective double subscripts, e.g., the gauge generator for the transformation synchronous- 
comoving to the total matter gauge is denoted with £g T . For an overview of our used notation 
see tab. 1 on page 33. 

2 The FLRW metric perturbed up to second order 

The components of a spatially flat FLRW metric perturbed up to second order are written as 


9oo = -a 2 (1 + 2^i + ^ 2 ) 

(2.1a) 

goi = a 2 ^Bn + - B2i^j 

(2.1b) 

Qij — & [&ij + 2 Cl ij + C*2 ij ] , 

(2.1c) 
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where the background part (which is by definition only time-dependent) is given by 

5oo = -a 2 0?), 9ij = a 2 (v) Sij ■ (2-2) 

Here, rj is the conformal time (a dr] = dt, where t is the cosmic time), and a(r]) is the FLRW 
scale factor, which obeys the Friedmann equation (see section 5.1). It is convenient to split 
the perturbations into the so-called scalar, vector and tensor parts, 


B (r) i ^(r) T ^(r) i (2-3) 

111 

C(r) ij 4 > (r )^ij T DijE (j.) + — diF^j -(- — 0jF] r ^ j -(- ~X(r)ij ) (2-4) 


where (r) = 1,2, stand for the rth order of the perturbation, we make use of the operator 
Dij = didj — (1/3 )'V 2 5ij, and W( r ) j and are transverse vectors, i.e. = d t Fr r \ i = 0 

and X(r)ij i s a transverse and trace-free tensor, i.e., d l X( r )ij = 0, x( r j i = 0. Historically, 
the reason why such a splitting was introduced [27, 43] is that at first order these different 
perturbation modes are decoupled from each other in the perturbed evolution equations, so 
that they can be studied separately. This property does not hold anymore beyond the linear 
regime, where higher-order perturbations are sourced by products of lower-order perturba¬ 
tions. Of course, despite the fact that different perturbation modes generally do not decouple 
at higher orders, the above splitting is still physically meaningful. 

So far our considerations are fairly general. In the case of the ACDM Universe, the 
metric in (2.1) can be simplified, as we can consistently neglect first-order vector and tensor 
perturbations: linear vector perturbations are not produced by standard inflationary scenarios 
and in any case decay with the expansion of the Universe [23], and tensor perturbations are 
believed to give a negligible contribution to the matter dynamics. However, we should keep 
in mind that the same reasoning does not apply to second-order perturbations: in the non¬ 
linear case, scalar, vector and tensor modes are coupled and the second-order vector and 
tensor contributions are generated by first-order scalar perturbations even in the absence of 
linear vector and tensor perturbations, [23]. With these simplifications, the spatially flat 
FLRW metric perturbed up to second order is 


9 00 = -a 2 (1 + 2^1 + V>2) 

901 = o 2 ( diB\ + — diE>2 + — UJ2 


9ij — ® 


(1 — 2</>i — (j>2) Sij + (2E\ + E 2 ) + -diF2j + —djF2i + —X 2 ij 


(2.5a) 

(2.5b) 


(2.5c) 


The four-velocity of matter is u^ = dx^/d r, where r is the proper (comoving) time, comoving 
with the fluid. For latter convenience we set 

„** = I ($ + „"), (2.6) 

where v^ /2 + ... is the peculiar velocity (peculiar in the spatial and temporal 

sense). From the normalisation condition = —1 we obtain the constraint for the 

time component of which reads up to second order (in any gauge) 

13 1 

v° = -Vq - 2^2 + 2^1 + 2 VlkV i + v i kB i ■ ( 2 - 7 ) 
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The perturbations of the spatial components v l split as usual in scalar and vector parts 


v \r)= S ' 3 V(r)j+w\ r) , ( 2 . 8 ) 

where d t w 1 ^ = 0. Finally, the perturbation in the matter density is written as p = p + pi + 
P 2/2 + ..where the background density p is time-dependent only. In this paper we work 
with the density contrast which is defined by 

(2.9) 

and is expanded as 

8 = 81 + y ( 2 . 10 ) 

up to second order. 


3 Gauge transformations in perturbation theory 


The theory of the gauge transformations in any given background space-time and beyond the 
linear order was applied to cosmology in refs. [32, 33], following the approach of refs. [39, 40]. 
See [44] for a recent review. 

We will adopt in the following the so-called passive approach, where a gauge transfor¬ 
mation is seen as a coordinate transformation —> x^. Up to second order it reads 


\x a ) = x» - g(x a ) + - Ci(x a ) - #(®°) 


with its inverse 


1 r 


x a {x») = x a + (#*) + - [^p(x») g (^) + & (s'*) 




(3.1) 


(3.2) 


where all the quantities are evaluated at the same point on the background space-time where 
the coordinates x^ and x 11 coincide. As usual, the four vectors can be decomposed into 
scalar and vector parts 


£(*r) = «(r), C( r ) = d l /3( r ) + d\ r) , with did\ r) = 0. (3.3) 

In order to find the components of a tensor in the new coordinates we simply start from the 
standard transformation rule for the metric tensor 


~ 9 ^ a ) = d^w 9 ° ^ Q( * P)) ’ ( 3 - 4 ) 

and expand the Jacobian matrix, and the argument of the metric components on the r.h.s. 
up to second order. The result of this infinitesimal coordinate transformation is naturally 
expressed by means of the Lie derivative, 2 C g. We have up to second order 


9/-US — 9fj.u + •AjS'/ii' T „ . 


"The explicit expressions for the Lie derivative can be found in appendix B. 
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(3.5) 









In general, by expanding the standard transformation rules under the coordinate transforma¬ 
tion (3.1) for any unexpanded quantity T - a scalar, a vector, or a tensor —, we find that 
its change up to second order is given by 

± = T + £ 6 T + ^(£| 1 T + £ 6 T) . (3.6) 

In order to find how the perturbations of a given tensor held transform, we consider now the 
expansion up to second-order of a generic tensor held T in two different gauges, 

T = T + Ti + ^T 2 , and T = T + Ti + ^T 2 . (3.7) 

The transformation rule for the perturbations is easily obtained by plugging (3.7) in (3.6) 
and collecting the terms of the same order. It reads 

T + Ti + ^T 2 = T + Ti + £ 6 T + ^ (T 2 + 2£ 5l Ti + C \ ± T + £ 6 T) . (3.8) 

This relation tells us how to obtain the perturbations in one gauge from the corresponding 
perturbations in the other gauge, that is 

Ti = Ti + £ Cl T (3.9) 

at hrst order, and 

t 2 = T 2 + 2%Ti + £| i T + £ 6 T (3.10) 

at second order. The r.h.s’s of these equations are automatically written in the new gauge. 

In the following two subsections, we give the transformation rules for the quantities we 
are interested in, namely the metric tensor, the density contrast and the four-velocity. They 
are straightforwardly obtained from eqs. (3.9) and (3.10) and by the use of the expressions in 
appendix B. 

3.1 First-order transformations 

Metric tensor We find the following first-order transformations for the metric tensor: 

• scalar perturbations 

Vh = Vh + 'Holi + ati 

0i = 0i — Hati — 

B\ = B\ — oi\ + Pi 

E\ = Ei + Pi, 

• vector perturbations 

&li = Wli + du 

Eh — En T du , 

• tensor perturbations 

Xiij = Xiij , (3-17) 

where the dot denotes partial derivative with respect to conformal time and 7~l = a/a = aH 
is the conformal Hubble parameter. Note that, since we neglect first-order vector and tensor 
perturbations, the first-order vector part of the spatial transformation d\ is constant and can 
be set to zero. 


(3.11) 

2 Pi (3.12) 

(3.13) 

(3.14) 


(3.15) 

(3.16) 
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Three-velocity The transformation of the temporal part of the peculiar velocity is obtained 
from eq. (2.7) and reads 

Vi = v® — T~Lai — a \. (3.18) 

For the scalar and vector part of the spatial peculiar velocity we find 

v\=v\- $1 (3.19) 

and 

w{ = w[- d\ . (3.20) 

Matter density Finally, the perturbation of the density contrast transforms as 

~5 l = 5 l -2,Ua l . (3.21) 

3.2 Second-order transformations 

In order to write in compact form the second-order gauge transformations we collect the 
contributions from products of first-order transformations defining II, Epflj and T ij. Their 
explicit expressions are given in Appendix C. 

Metric tensor The second-order transformations for the metric are: 

• scalar perturbations 

fa = fa + Ha.2 + 02 + II (3.22) 

fa = Ua.2 ~ lj k k ~ \v 2 fa (3.23) 

fa 3 

B-2 = B 2 - 02 + fa + V~ 2 S fc fc (3.24) 

E 2 = E 2 +fa + ijv- 2 v- 2 r^. - ^V- 2 T\., (3.25) 

• vector perturbations 

fa i = ^2 i + fai + Sj — V (3.26) 

F 2 i = F 2 i + d 2i + V _2 Tj fc fc - V- 2 V~ 2 r k \ kH , (3.27) 

• tensor perturbations 

fan = X2ij + faj + - 2 (v~ 2 r k \ kl - t \.) Sn + lv- 2 v~ 2 r k \ kHj + 

+ \v- 2 T k Kij - V- 2 (X 4 %. + , (3.28) 


where V 2 stands for the inverse of the Laplacian operator with Euclidean metric. 

Four-velocity From the normalisation condition (2.7) in the new gauge we find the trans¬ 
formation of the temporal part of the peculiar velocity which reads 



v 2 = v 2 ~ F.O-2 — fa + Oi\ 2 (yi — Hv °) + ai (jd . 2 — 7 'ij + Tiai — di 



+Ci (, 2v i ,i ~ - dp,:) + di (&i - 2v°) - 2a 1}i v\ + aifa\ . 

(3.29) 

The scalar and vector part of the spatial peculiar velocity are given by 



V2=V2-fa + V“ 2 H fc fc 

(3.30) 

and 

w i 2 = w i 2 -di + n i -v~ 2 n k M . 

(3.31) 
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Matter density Finally, the second-order perturbation of the density contrast transforms 
as 


&2 — 3’H«2 + CX 1 


6W 


+ 2 


9 Ufa 


'm 0 


— 31~Lai + 25i — QH5i 


+ (25 1 -3Ua 1 ) k ^. 

(3.32) 


4 Lagrangian gauge and Eulerian gauges 

The Lagrangian frame for irrotational and pressure-less matter, in the language of gauges, can 
be associated with the synchronous and comoving gauge (see, e.g., [1, 6, 7]). The synchronous 
gauge is defined by [42] 

Bg = = 0, (spatial gauge condition) (4.1) 

ipS = 0. (temporal gauge condition) (4-2) 

In such a coordinate system, the fluid particle is at rest which amounts to use the comoving 
gauge conditions vq = Wq = 0. 3 We note that the simultaneous conditions “synchronous” 
and “comoving” hold only for an irrotational and pressure-less fluid, and when neglecting any 
residual gauge modes [28, 44, 46]. 

From the quote in the introduction one can deduce a general recipe to find any Eulerian 
gauge. A major subset of Eulerian gauges are defined with the spatial gauge condition [6] 

E = F l = 0. (4.3) 

One reason for chosing this spatial gauge condition is the following. In the Lagrangian frame, 
the dynamical information of the displacement field is encoded in Eg and Eg. Therefore, 
when performing a spatial gauge transformation to an Eulerian frame with vanishing E and 
F l (see eq. 4.3), the spatial gauge generator must carry this dynamical information. In that 
very case, the spatial gauge generator is the displacement field. 

Having specified only the spatial gauge condition is of course not sufficient for the study 
of general relativistic perturbations. A temporal gauge condition has to be imposed, which 
is, in principle arbitrary: Any temporal gauge condition supplemented with the spatial gauge 
condition (4.3) fixes a unique Eulerian gauge. Then, the resulting Eulerian equations of 
motion take, at least to first order, the Newtonian form in Eulerian coordinates [19]. 

In this paper we shall derive Eulerian solutions for two gauge choices, namely for the 
Poisson gauge and the total matter gauge. The Poisson gauge is defined by 

Ep = Fp = 0, (spatial gauge condition) (4.4) 

.Bp = 0. (temporal gauge condition) (4-5) 

It was originally introduced in ref. [28] as representing a generalisation of the so-called New¬ 
tonian (or longitudinal gauge) in order to account for vector and tensor modes, which are for 

example generated by the coupling of scalar modes beyond the linear regime. Note that the 
temporal condition corresponds to the minimal shear hypersurface condition in ref. [43] (see 
also [28]). 

3 In this paper, we do not explicitly use the comoving gauge condition but rather use it as a consistency 
check for the derivations in section 6. Nevertheless, also since this is common in the literature, we shall call 
the gauge corresponding to the conditions (4.1)—(4.2) the “synchronous-comoving gauge”. 
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The total matter gauge is another Eulerian gauge which is Eulerian in terms of the spatial 
coordinates, i.e., the spatial gauge conditions are given by E7p = Ep = 0. The temporal gauge 
condition is defined by requiring the vanishing of the total momentum potential on spatial 
hypersurfaces (see, e.g., [44, 45]), i.e., the vanishing of the scalar part of the °j component 
of the stress-energy tensor T IJ U . At first order, the perturbation of the total momentum is 
T\i = {pF P) (vu + Bu), where P is the background pressure, and the resulting temporal 
gauge condition reads 

ui T + Bi t = 0 . (4.6) 

At second order, the total momentum is [44] 


T^={p + P) 

+ 2 (pi + Pi) (vu + Bu) + i^B\ + Vi^j link j 


V2 i + B2i + / iCukV 1 — 2-01 (vu + 2Bu) 
2 
a 2 


(4.7) 


where II ij is the anisotropic stress tensor. Assuming vanishing first-order vector modes and 
the spatial gauge condition of the total matter gauge, the requirement of the vanishing of the 
scalar part gives at second order 


V2 t + B2 t — 2V 2 d l (Vh T E>i T ,z + 2(j>i T vi Tt i) = 0 . 


(4.8) 


This temporal gauge condition is valid for a (multi-component) fluid with pressure and non¬ 
vanishing anisotropic stress. Furthermore, we note that there exists an alternative definition 
for the temporal gauge condition of the total matter gauge [24, 25, 46], that is employing 
the vanishing of the scalar part of the spatial component of the covariant 4-velocity, m = 
g i0 u° + gijVp . Demanding m = 0 for its scalar part leads to the same condition (4.6) at first 
order, and to condition (4.8) at second order, 4 and thus this temporal gauge condition is 
formally equivalent with demanding the vanishing of the scalar part of T°j. In summary, we 
define the total matter gauge with 


Ep = F?p = 0, (spatial gauge condition) (4.9) 

S l T°i T = 0, (temporal gauge condition) (4.10) 

where we have defined the scalar mode extraction operator 5* := V _2 <9* and, again, condi¬ 
tion (4.10) amounts to use (4.6) at first order and (4.8) at second order. 5 As we shall see in 
section 7 when we calculate the perturbations in the total matter gauge for an irrotational 
dust fluid, the time coordinate in the total matter gauge is the proper time of the fluid parti¬ 
cle, as it is the case in the Lagrangian gauge. In the language of gauge transformations, this 
also means that the temporal gauge generator from synchronous-comoving gauge to the total 
matter gauge is precisely zero. This choice of Eulerian frame is thus very convenient when 
relating GR results to Newtonian investigations, since the relativistic Eulerian-Lagrangian 
correspondence makes use of the identical time coordinate, i.e., the proper time. So to say, 
the problem of relating a GR description to a Newtonian world is reduced from a 4D problem 
to a 3D problem. Such considerations are for example useful when generating GR initial 
conditions for Newtonian IV-body simulations (see [47]). 

4 We thank John Wainwright and Claes Uggla for pointing this out. 

5 In an earlier version of this manuscript, we have used the temporal gauge condition v° = 0 instead 
of (4.10). We stress that for an irrotational and pressure-less fluid the results in the total matter gauge are 
unaffected by choosing that condition or the other, yet it is highly non-trivial to prove the equivalence of these 
conditions on general grounds. We leave this open issue for future work. 
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5 Evolution of perturbations in the Poisson gauge 


The Poisson gauge can be perturbatively defined by -B( r ) p = 0 for the temporal condition and 
by E ( r \ p = -E(r) P = 0 for the spatial condition. For second-order perturbations in the Poisson 
gauge in the Einstein-de Sitter (EdS) background see e.g. [23] and [29], and for a non-vanishing 
cosmological constant see e.g. [22, 30]. Here we start by reporting and simplifying the results 
of ref. [22] for the metric and the four-velocity, and of ref. [11] for the density contrast. The 
new expressions obtained here will be the starting point for the gauge transformation in the 
following sections. Let us just recall the procedure to solve the Einstein equations for the 
perturbations, which is the same at first and second order, see Appendix A of [22], The trace 
part of the ij components of Einstein equations gives the evolution of the scalar perturbation 
0( r )P, and the trace-free part gives the relation with the other scalar ^(rjP- The vector mode 
0 J 2 i P is found from the 0 i components of Einstein equations and the tensor modes 7T2 ij P obey 
an evolution equation obtained from the trace-less ij components, once the scalar and vector 
perturbations are determined. Finally, the 00 component of the Einstein equations and the Oi 
components, together with the momentum conservation, provide the density and four-velocity 
perturbations, respectively. 

Following [22], we write the initial conditions in terms of the curvature perturbation of 
the uniform density hypersurfaces, C- This is gauge invariant quantity and remains constant 
on super-horizon scales after it has been generated at the primordial epoch [26], and therefore 
provides the initial conditions for cosmological perturbations re-entering the Hubble radius 
at later time, including all the necessary information about primordial non-Gaussianity. We 
expand ( as C = C^ + (l/2) £( 2 ) = C^ + (a n i — 1)(£^) 2 + - ■ •, where the parameter a n \ encodes 
the local primordial non-Gaussianity for different inflationary scenarios. Initial conditions are 
fixed deep in the matter-dominated era when the relation between the curvature perturbation 
with the gravitational potential is given by Cl = — 5<£>ii n /3. 6 Our initial conditions up to 
second order are then written as 

1 5 25 

Cin = Clin T 2 < ’ 2in = — ^ g~ (® n l _ ^ VTin i (5*1) 

where Clin and ipu-n. are first-order Gaussian random fields. The same initial conditions are 
written in terms of the initial gravitational potential as ip i n = ipn n + f n \ pf in , where we can 
easily read-off the relation of the non-linear parameters, 

fni = ^ (a„i - 1) . (5.2) 


5.1 FLRW background equations 

Let us briefly summarise the results of the evolution equations for the FLRW background. The 
energy constraint, Raychaudhuri equation and continuity equation for matter are respectively 


3V 2 = SirGa 2 p + a 2 A , 

(5.3) 

3V. + AirGp — a 2 A = 0 , 

(5.4) 

p + 3 Vp = 0 , 

(5.5) 

°This relation is sometimes also given in terms of the Bardeen potential 
gravitational potential is $ = —<p lin . 

Its relation to the linear 
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where A is the cosmological constant and p ~ a 3 is the background density of matter. We 
will make frequent use of the following definitions 


_ SttGu 2 Ultimo 0 _ a 2 A _ o 

* m 3% 2 P aH 2 ’ A_ 3 K 2 U 2 ’ (5 ' 6) 

where the subscript 0 indicates an evaluation at present time r/o, and we have defined fl m o = 
87rGagPQ/(3'Ho). By using these definitions the Friedmann equations (5.3) and (5.4) read 

V 2 = Ul (fw 1 + fW) , (5.7) 

n = n 2 -l^h^. (5.8) 

2 a 

In this paper we will always assume a flat Universe with O m + Oa = 1. 


5.2 First-order results 

Metric tensor Since we assume vanishing vector and tensor perturbations at first order, in 
this section we only consider scalar perturbations. At first order, the trace-free part of the ij 
components of the Einstein equations implies that the two scalar perturbations in the metric 
coincide, i.e., 

Vfip = 0i P = <P, (5-9) 

and the trace part of the ij component gives the evolution equation 

(p + ‘i'Hj) + a 2 A<p = 0 . (5.10) 


The solution, selecting only the growing mode, can be written as 


V(x,rj) = g(rj) (p 0 (x), (5.11) 

where ipo is the peculiar gravitational potential linearly extrapolated to the present time r/o 7 
and the time coefficient g has still to be determined. Plugging the solution (5.11) in the 
evolution equation (5.10) and introducing the new unknown variable V = ag, we obtain the 
differential equation 

V + llV- _ ^o° m0 £> = o, (5.12) 

2 a 

which is nothing but the known evolution equation for the first-order time coefficient of the 
Newtonian density contrast [34], In this paper we are only interested in the growing mode 
solutions of V (and of the other time coefficients, of course). The solution for the first-order 
scalars is thus given by (p(x, rj) = g(rj) </?o(*), where g = T>/a is identified to be the growth- 
suppression factor, and V is the growth factor, namely the linear growing-mode solution of 
eq. (5.12) for the Newtonian density contrast 5i(rj,x) = T>(r])di(x,go). We note that this 
identification implies that the solution in eq. (5.11) coincides with the first-order Newtonian 
cosmological potential. The analytical expression for the fastest growing mode V in terms of 
the scale factor is given by [35, 36] 


, 5„ (2 „ H f a dx 

V(a) = -U^ m0 - 


— a\ 1 + 


^A0 


3 5 11 U A0 


a 3 2 A 1 

1 2 6 6 




(5.13) 


1 Here and in what follows quantities with the subscript 0 are meant to be evaluated at the present time. 
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where 2 -Pi is the Gauss hypergeometric function. We normalise the growth factor such that 
T>(r] 0 ) = 1. We finally consider the Newtonian Poisson equation 


V (p(x,rj) - - 


3 nln 


0 s ^mO 


Si(x,rj) = 0, 


2 a 

evaluated at present time to fix the constant d and we find 

jf , ^ 2 V 2 <po(x) 

h(x,Vo) = 0^727, -• 

o rtgiimO 

In an EdS universe we have / D(rj) = rf and = 4, so = (r/ 2 /6) V 2 </?o- 

Four-velocity The spatial component of the peculiar four-velocity is given by 


(5.14) 


(5.15) 


1 


vip = - 


n (V + Ucp) = -- 


47r Ga 2 p x ^ ' ' 39-Lln 


-Vipo , 


0 12 m0 


(5.16) 


which contains the scalar part only (the vector part of the velocity is vanishing at this or¬ 
der). This coincides with the first-order Newtonian peculiar velocity in Newtonian Eulerian 
perturbation theory (NEPT). To obtain the perturbation of the zeroth component of the 
four-velocity, we plug the result (5.11) in Eq. (2.7) and immediately obtain 


^i P = -Vh P = ~g <Po ■ 

Matter density The final expression for the density contrast is 

2 


£i P — 


min 


0 14 m0 


wipo — mvipo 


(5.17) 


(5.18) 


The first term in the above equation is identical to the first-order Newtonian density contrast, 
whereas the second term is a GR correction. Actually, as we shall see in section 6.1, this GR 
correction is the result of the different time coordinate used in the Poisson gauge, which is 
due to the relative velocity between this Eulerian frame and the Lagrangian frame, see in 
particular eq. (3.21) and eq. (6.3). 

5.3 Second-order results 

Metric tensor: scalar perturbations Our starting point are the results of [22] and [11], 
However, as we shall see, it is possible to simplify their expressions tremendously. We begin 
with eqs. (2.22) and (2.23) of [22] 


10 


■ 02 p ( 77 , x) = (Bi+ 4 g 2 - 2gg in - — (o n 1 - l)ggin ) Vo + 6 




3 

e+ 2 


•,99'va 


00 + £> 3 V 2 d i d :, (d l (podjipo ) + , 


10 


4> 2p(v, x ) = [Bi - 2 99in - y ( fl ni - 1 )99mj + 6 ( B 2 - -ggin ) ©0 

+ B 3 V~ 2 did J {d l (f 0 djip 0 ) + B^tpodiipo , 


(5.19) 


(5.20) 
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where the subscript “in” denotes an evaluation at initial time (given deep in the matter 
dominated era), e = / 2 /n m with / = dlnP/(dlna), and 



1 

3 


Vo <Po,i 


7-2 


,/ ,m 

Vo Vo 


,/m 


(5.21) 


The time coefficients Bj(r]) are reported in [22] in eqs. (2.24) and (2.25) (in the notation of [22], 
these coefficients are Bi). They involve unfortunately very complicated integrals over time 
with no explicit solutions given. Here we find explicit solutions for these time coefficients. Let 
us begin with the PDE for the scalar perturbation <^> 2 P which is the trace of the ij components 
of the Einstein equations, eq. (2.7) in [22] 

<j> 2 p + 37~L<j)2p + a 2 A</> 2 p = S(rj,x ), (5.22) 

with the source term 


S(rj,x ) = g 2 Vi ra H 2 


(/- l) s 


+g 


n D 

^ (jr + (djVo&Vo) - ditpo&’ipo 


Vo + 12 2 


(/-l ) 2 


~st + 3 00 


To obtain the ODE’s for the Bi s we impose the Ansatz 
4>2p iv, x) = B\{ri) <(4(*) + B 2 {rj) 4p(*) + ^(v) 4>[p( x ) + 4? (*) + — 

Hin 


02p 


(5.23) 


(5.24) 


(5.25) 


where 

4p (*) = <Po » 4? (*) = ©0 , 

4p (*) = V~ 2 9j(9 J (d'tpodjtpo ), 4p 0*0 = d l VodiVo , 

and the last term is the solution of the homogeneous equation with d> 2 (x) representing the 

initial condition taken deep in matter-dominated era and on super-horizon scales. We now 
repeat the same procedure we used to determine the time coefficient of the first-order scalar 
(j>i p , given in eq. (5.11). For demonstrative purposes, we briefly show how to solve for the B\. 
We plug Bi<f> 2 p into (5.22) and keep only terms proportional to 4j? = <pQ. The temporal and 
spatial dependence thus factorise so that the ODE for B\ is 

Bx + 3nBi + a 2 A£q = g 2/ U 2 {f - l) 2 , 


or, introducing the new unknown variable b± = aB±, equivalently 


bi+Ubi-- 


3 nln 




a a 

where we have used eq. (5.8). Likewisely, we obtain for the other bi = aB, L 
b 2 + m 2 - l^h^b 2 = ^44 [2(/ - l) 2 + 3(H m - 1)] 


a a 

v 3 ^fi m0 , V 2 

03 + 7ib 3 - -- b 3 = 2-b 

2 a a 

i i -l ju 3 77gH m o^ _ T>' 2 

2 a a 


4T> 2 


min 


0 s ^mO 


(5.26) 

(5.27) 

(5.28) 

(5.29) 

(5.30) 
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To arrive at eq. (5.29) we have used the ODE for the first-order growth D, eq. (5.12). The 
easiest solution is obtained for 64 , by noting that its ODE can be written as 


j: \ 'Ylj t 3'Hg^niO j- _ 3 475D m 0p2 

2 0 2 a 

(5.31) 

with T = — (3/2)%gfi m o&4- Equation (5.31) is nothing but the known evolution equation for 
the second-order time coefficient of the Newtonian displacement field . 8 Its solution for the 
fastest growing mode in a ACDM Universe can be found in [37]. For an EdS universe we have 
T = (3/7)r/ 4 . The solutions for 61 , 62 , and 63 are also easily obtained. We find the fairly 
simple growing mode solutions 

, ^ 2D 2 

61 - Va+ ■3Hin m „ + 3 Vsh "' 

b 2 = -2V(g in - g ), 

(5.32) 

(5.33) 

63 - 3Wgn m „T + v ) • 

(5.34) 

64 =- 3 - J- . 

3%oD m o 

(5.35) 


Note that for an EdS universe, b\ = 62 = 0, whereas 63 = (5/21)r ^ 4 and 64 = — (l/14)r/ 4 . 
With the above simplifications, we can finally write the two scalar perturbations as 


V-% = 3 g 2 + -c/5in(l ~ 2a ni ) + 


2D 2 

,1 4(D 2 +D) 


, „ ( . 2 10 4 V 2 

*’ 5 + 6 49 - J Sgi “ + 3 M ^ 


2D 2 

+ 3a^D m0 lf ° m ~ 3 aUltt 


Tr 


4>2 p = \-g 2 + -P5in(l - 2a n i) + 


2D 


52 


0 1 ‘mO 
2 


+ 


2D 2 


3a?^QD m o 


3a%gD m o 
,/ 4(P 2 + P) 


10 


^0 + 6 2 £f - —fif0i n @0 


¥> 0,/¥> 0 


Tr 



(5.37) 


where 


d'n = -V- 2 - 


(V (p 0 ) - <po,ik<Po 


,ik 


(5.38) 


The last line in (5.36) and (5.37) is the second-order gravitational potential in NEPT (see, 
e.g., [37]), whereas the remnant terms are relativistic corrections. 

Note that in deriving the last line of (5.36) (and the last line of (5.37)), we have made 
use of the identity 

V~ 2 did rn {d l pod rn ipo) = -2T 0 + • (5-39) 

8 Note explicitly, that J- is not the second-order growth of the density because the density does not factorise 
into a single time and space-dependent part. By contrast, the displacement field does: In the notation of 
eq. (1.1), the Newtonian displacement field is up to second order 4q = u(q) + T 2 i(q), where V and 

T are truly the solutions of eqs. (5.12) and (5.31), respectively. For the Newtonian trajectory see eq. (8.19). 
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Using this identity in (5.37) we obtain the relation between the “GR kernel 11 ©o and the 
“Newtonian kernel” \ko: 

V 2 ©o = d'o - 2 ^o,/¥ 3 o > (5.40) 

which will be helpful in the following sections. 

For the limiting case of an EdS universe with a n i = 0, our results (5.36) and (5.37) agree 
with [23]. For ACDM, the expressions (5.36) and (5.37) have been first derived in refs. [24, 25]: 
Although different spatial and temporal functions 9 are used, our result agrees with theirs. 


Metric tensor: vector perturbations In a ACDM Universe, linear vector perturbations 
are usually discarded. But even in the absence of linear vectors, the non-linear gravitational 
evolution generates second-order vector perturbations, as it is well-known. These second-order 
vector perturbations arise from quadratic couplings of first-order scalar perturbations. For an 
EdS universe, the first derivation of these vector perturbations was given in [23], whereas the 
results for a ACDM Universe were reported in [22], To get these vector perturbations, one 
possible derivation is to first observe that the second-order momentum conservation equation 
= 0. see e.g., eq. (A.26) in [22]) is only sourced by a pure scalar for growing mode initial 
conditions. From that one can easily deduce that the vector perturbation uj 2 i P has to be equal 
minus the transverse part of the velocity, —W 2 i P - Then, UJ 2 i P can be extracted from the 07 
component of the Einstein equations which is [22] 10 


1 


77-0 2 p,i + 02p ,i ~ -V W 2 ip + 10gg(p 0 ipo t i = -377 D m (gtp 0 + <5i p )-ui p ,j + 


V2 P ,i 


(5.41) 


Extracting from this expression the transverse part we then obtain 

16 VV 


W2i P — ~- 


-77 i(x) 


377oD m o a 

which is purely relativistic, and we have defined the transverse kernel (V ■ 77 = 0) 

TZi(x) = V -2 (<A),;VVo - ¥>bVo,Zi + 2^ 0ji ) . 


(5.42) 


(5.43) 


Metric tensor: tensor perturbations In the Poisson gauge, the tensor perturbations 
obey the wave equation [22] 

K 2 ij p + 2777T2ijp - V 2 ir 2 ij p = s(rj) %j(x) , (5.44) 

where s(rj) = —8g 2 (l + 2/ 2 /(3f2 m )), and Tij(x) = V _2 5jj is a divergenceless and traceless 
tensor related to 


Sij{x) — V 'kg dij + ~b 2((^o, ^POjik^P’oj) ■ (5.45) 


The solution of (5.44) is 


/ d 3 fc 

j^TT2ij P (ri,k)e lk ' x , (5.46) 

9 In particular, their temporal function B is related to our T and T> via B = (1/3)[1 — 2 TI'D 1 ]. We thank 
John Wainwright and Claes Uggla for pointing this out. 

1() Note that in deriving this expression we have corrected two typos in eq. (A.27) of [22]. 
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with 


(5.47) 


^2ij P (v, k) = %j(k ) h(r], k) 


where k = \k\, and 


n 


Mv,k) = Xi(jj,k) [ dfi^%^-s{fj) ~X2(v,k) [ d 


W 

where xi and X 2 are the solutions of 

X + 2Hx + k 2 x = 0 , 

with corresponding Wronskian W = X 1 X 2 ~ XiX.2- 


W 


(5.48) 


(5.49) 


Four-velocity To our knowledge there are no explicit solutions in the literature for the 
second-order velocity in the Poisson gauge for the ACDM Universe. We take the divergence 
and the curl from eq. (5.41), then we obtain respectively 


2 a 


V2 P = ~ 


and 


3?7gU m 0 
16 VV 


(h^2 p + </> 2 P + 5gg<Po) 

*0 + 


2 gV _ 2 mV 1 
+ 3^U m o^ 0 “ 3 {-Uprao) 


2^0 


9(^^mo) 2 ^ + 9(H 2 0 n m o) 2 ^ 0 ^ 


4VV 


2V 


3^U m0 L 

4 VV 


10 

T 


9 + ^ffin(Onl - 2) ipl 


,1 


9 (^Ono) 2<w ° + 9 (n 2 n m0 y 


T 


'W2iv> ^2zp 


^OnO 00 


d'o, 

(5.50) 

VV 

(5.51) 

Hi, 

> ^ L ’ 


where we remind the reader that the transverse kernel lZi(x) is defined in eq. (5.43). Note 
that the last line in eq. (5.50) is the second-order peculiar velocity in NEPT (see, e.g., [37]). 
To arrive at the final expression for V 2 P , we have made use of the identity 


nv 


5 

— n^in 


n 2 n m0 2— 2 


9- 


(5.52) 


which we shall prove later, see in particular the derivation of eq. (6.13). For the perturbation 
of the zeroth component of the four-velocity, we plug the result (5.36) in Eq. (2.7) and obtain 


U2p = - I ^99i 11 (1 - 2o.nl) + 


2V 2 


+ 


4V 


52 


3a77oU m0 

2V 2 ^ 


2 / , , 10 4X> 2 

v ° - 6 1 ■ 49 - T 99 " + 


©0 


9 {H 2 n m0 ) 2 3 an 2 n m0 


1 4 (V 2 + F) 


(5.53) 


We have checked that the EdS limit of the above results agree with the ones in [23, 31]. 


1 Note that we have corrected a typo in the last integrand of eq. (4.11) in [22]. 
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Matter density With the explicit solutions of the time coefficients fej given in eqs. (5.32)- 
(5.35) we are able to simplify the result given in ref. [11] 12 to the fairly compact expression 


2'h 2 t > 2 r „ 2 2 , ion „ , 0 . ■ 2 2mvv 

!p = ^ " 4/ ] + (1 + 2 a ni)^in^o - nV 2 n Qo 

art-O^mO on.Qi£ m o OrTQiimO 


+ 


+ 


+ 


2 £> 20 ^ 
n 9 7 t9 inO n l 

^Q^rnO L 9 


>2 87-// 

(Vw) + 3(«Sn m0 )^° 


37/2P 


■O^mO 

4 


10 4 

6 Pg + —Vg in {\ - 2a n i) + - 


V 2 


3 7/ 2 54 


O^mO 


<A)VV 0 


9(7/ 2 fl m n) 2 


(D 2 + ^)(V 2 ^o) 2 + 2P 2 ^V 2 (, 




(5.54) 


We have checked that the EdS limit of the above agrees with the results of refs. [23, 31]. For 
ACDM, our result is in agreement with the ones reported in [24, 25], however care has to be 
taken when comparing our results to theirs, because of different notations and normalisations 
used; for a clarification of this we refer the reader to section 4.4.3 in [25]. Note that the last 
line in eq. (5.54) is the second-order density contrast in NEPT (see, e.g., [37]). 

Finally let us conclude this section with a comment about our solution scheme for 
Eq. (5.22) and its applicability to cosmological models beyond ACDM. The general solu¬ 
tion technique for cj> 2 P (which comes with the Ansatz (5.24)) should also be fruitful for other 
cosmological models, as long as the source term S of the PDE for (j) 2 P , Eq. (5.22), can be 
written in factorisable form. Note however that the source term cannot be factorised in e.g. 
multi-fluid models, or when a radiation component is included in the analysis, because then 
the respective growth functions are generally scale-dependent (see e.g. [38] for the massive 
neutrinos case). 


6 From the Poisson gauge to the synchronous-comoving gauge 

In this section we apply the general formulas of section 3 to the transformation from the 
Poisson gauge to the synchronous-comoving gauge, which is defined by il>( r ) s = 0 for the 
temporal condition and by -B( r ) s = w i(r) s = 0 f° r the spatial condition. We define the 
synchronous-comoving gauge as the Lagrangian frame of reference of the fluid particle. The 
respective Lagrangian coordinates are denoted with (r, q). 

Let us note that the equations for the metric transformed from any gauge to the 
synchronous-comoving gauge have always the following structure: at any fixed order the 
equations for the scalar in goo an d f° r the scalar and the vector in goi are f° ur equations 
for the unknowns «( r ), /3( r ) and d 1 ^. i.e., they give the functions of the transformation . 
These results then give the spatial metric in the synchronous-comoving gauge. The equations 
for a( r ), /3( r ) and dl, contain time derivatives, and the integration constants in the results 
represent the residual gauge freedom of the synchronous-comoving gauge. It can be set to zero 
by requiring that at initial time spatial and temporal coordinates coincide in the Lagrangian 
and Eulerian gauges, i.e., the Lagrangian and Eulerian positions coincide initially. We shall 
make frequent use of this requirement. 

12 Note that the first and the last coefficients of <^ 2 = g/o in the first line of ref. [11] eq. (29) are wrong. 
They have to be replaced by2(/ — l) 2 . 
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Before beginning with the calculations, let us emphasise that our results for the gauge 
generators and for the metric in the synchronous-comoving gauge are entirely new and they do 
agree with the findings in the literature in the restricted case of an EdS universe, [23, 41]. In 
addition, by transforming our result for the density contrast in the Poisson gauge, eq. (5.54), 
we are able to derive a fairly simple result for the density contrast in the synchronous-comoving 
gauge. Note especially that we are able to solve for all time coefficients which actually coincide 
with the first-order and second-order structure growth of the Newtonian displacement held, 
respectively denoted with V and T (see footnote 8 on page 14). 


6.1 First-order transformations 


Metric tensor From eqs. (3.11)-(3.14) we find 

0 = ip i p + Tloiip S + di PS , 

1 2 

01 S 0ip T~icx Ipg - V /3 Ipg , 

dips — /^lps > 

£ls — /^lps • 

From (6.1) and (6.3) we find the solution for ai PS and /?i PS 



( 6 . 1 ) 

( 6 . 2 ) 

(6.3) 

(6.4) 

(6.5) 

( 6 . 6 ) 


where the space-dependent integration constants I\ and I 2 represent the well-known residual 
gauge freedom of the synchronous-comoving gauge which is fixed from the initial conditions, by 
demanding that spatial and temporal coordinates coincide at initial time. Thus I\ = I 2 = 0. 
Using the identity 

^ 3M|SW-> 

2 « U n 

which is just the first integral of the equation (5.12) for the growth factor, we get 


a ips 


2 Vipo 


3 nisi 


s m 0 

o _ 2 *>¥*> 

P 1 PS — o n /2( 


3 H 2 0 n m0 

Finally, we find (pi s and E\ s by substituting the last two expressions in (6.2) and (6.4) 

A _ Vipo , 

01s 0-1/20 3" 9^P0 Ipg 1 

o rt-Qi imO 
rp 2 ^VO 

Elg - --- 


3 urn 


0 s 4 mO 


In our expression for <^>i s the quantity 


gipo - "Haips = \ g + 


2 nv 

S'Ho^mO 


<^0 


( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 
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is constant, as can be easily verified by taking the time derivative, and then using eq. (5.12) 
and the background Raychaudhuri equation % — T~L 2 + (3/2 )7-L 2 Q m = 0. That constant can 
be obtained by evaluating eq. (6.12) at initial time. We have thus 


<~p T^Ctlps ^Pin i 
where p- m = g-mPo- Then, our final expression for cf> i s is 


<£is 


2PVVo 5 

9 H 2 fl m o + 3 9inLP0 


(6.13) 


(6.14) 


Finally, by putting the above results together, we obtain the first-order result of the 
spatial metric 

,2 


9ijs ~ a 


10 \ 4 V 

3 PinJSij 3^2n m0 ^« 


(6.15) 


Four-velocity The transformation of the four-velocity gives simply Uj* = (l/a)(l, 0), which 
is expected since the peculiar velocity (both temporal and spatial parts) vanishes by definition 
in the synchronous-comoving gauge. 


Matter density Finally, the transformation of the matter density perturbation in eq. (3.21) 
gives for the density contrast 


_ 2 PV 2 y> 0 
ls “ 3 HfomO ’ 


(6.16) 


i.e., the first-order density contrast in the synchronous-comoving gauge is the solution of the 
Newtonian equation 


2 a 


(6.17) 


as it is well-known. The first-order relativistic correction in the Poisson gauge, which arises 
because of the non-vanishing velocity, is eliminated by the time transformation to the comov¬ 
ing (i.e., Lagrangian) frame. 

All of our first-order results agree with those of refs. [11, 26] for ACDM. 


6.2 Second-order transformations 

Metric tensor: scalar perturbations From the transformation rules at second order 
(3.22)-(3.25) we find 


0 — V^ 2 p T T~Loi 2 PS + 02 P g + Ilpg , (6.18) 

4>2 s = </>2p - T-La.2 p S - -T k kp3 - -V 2 /3 2ps , (6.19) 

0 = — 02 PS + $2 PS + V 2 T, k kps , (6.20) 

E 2s = /3 2ps + ^V- 2 V- 2 T^ .. ps - ^V- 2 T\. ps , (6.21) 


where the explicit expressions for the first-order squared terms n ps , Sj ps and T.;j PS can 
be found in Appendix C, whereas '(/ 2p and ^> 2p are given in eqs. (5.36) and (5.37). From 
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eqs. (6.18) and ( 6 . 20 ), it is then straightforward to obtain the following expressions for the 
gauge generators 


a 2 ps 


& 2 ps 


i n 


a 

rv 


dfj a 


Vin 


t rip 


dr) 


'Vin 


ct 2ps - V“ 2 S fc 


kps 


( 6 . 22 ) 

(6.23) 


To keep track of the computational steps, let us split a 2 PS into a “Newtonian”-like and “GR”- 
like part, i.e., a 2 PS = H ere i with “Newtonian”-like contributions, we mean that 

is identical with the Lagrangian velocity potential . 13 Indeed for the Newtonian part we 
have after some manipulations 




1 4T 0 

a 377/rimO 


dr) ( V 2 + V) 


'Vin 


9 {'H 2 0 n m0 ) 2 ’ 


(6.24) 


where the last expression is found by rewriting the differential equation (5.31) for the second- 
order growth function T in integral form: 


j- _ SjtfflrnO l\ n{p 2 +J:) 
2 a 


'Vin 


For the GR part of a 2 PS we obtain after some straightforward calculations 

1 n 


n- GR - 
a 2 PS - 


dr) 


' Vin 


1(1 . o o 

-7rT^9va\fl nl 1)^5 8 2 

o /l 0 li m 0 


/ 2 -/+^ n 


©0 


(6.25) 


(6.26) 


To arrive at the terms proportional to the spatial kernel ©o we have used the relations (5.6) 
and (5.52), and the fact that V = f'H'D. The integral proportional to is trivially solved 
by using (6.7). To solve the other part of the integral (6.26), we note that 


VV 

a 


U> w 2 p 2 ( / 2 - /+ N' 


which can be easily proven by applying the “reverted” partial integration rule 


1 

a l 


VV 


v 

Vin 


rv 

Vin 


dfj ( 


V 2 + vv 


(6.27) 


(6.28) 


discarding a decaying mode ~ 1/a, and by using the differential equation (5.12) to get an 
expression for V. Putting the Newtonian and GR part of the temporal gauge generator 
together, we then obtain 


2V 


CK2 PS 


377/Q m () 


y 5in (a n l - 1) vl ~ 12(70o 


+ 


S-TTo 


9 {n 2 n m0 y 


(6.29) 


Now we proceed with the calculation of the spatial gauge generator (6.23), which we again 
split into a Newtonian and GR part, / 32 PS = P^ps + ^2ps‘ For th e Newtonian part, the 
integration in time is trivial so we obtain immediately 


= 


9 (n 2 n m0 y 


- -V 2 (V (/? 0 ) 2 


13 For the Newtonian non-perturbative approach to this gauge transformation see [17]. 


(6.30) 
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where the last term originates from the Newtonian part of Ej ps , given in (C.6). For the GR 
part of /32 ps we obtain 


oGR = 

^2 PS Q-l/2f 


dfj 


Vgin ( ^a n i - | ) Vo - Vgipl - 6VgQ 0 


3 ^ 0^ m 0 JVin 

The time integral of the first term is trivial, and to solve for the latter we observe that 

r V 2 5 

/ d rfDg = 2 + -Vg in , 

Jrji n 6?^q fl m o 6 


(6.31) 


(6.32) 


which can be obtained from the identity T> 2 = 2 f d fj'D'D and then using (5.12) and (5.52). 
We then obtain the GR part of the spatial gauge generator 


« = 0172 




52 


10P<?in (flail 2) ^ /2 


2> 




V>0 


3’HnflmO 


P 


52 


77nfl m O 


+ 5Pg ir 


©o, (6.33) 


where we have set an integration constant to zero because we require that the Eulerian and 
Lagrangian frame coincide at initial time. 

We have thus solved for /32 PS = /S^g + {$2^- Equipped with the (longitudinal) gauge 
generators, we are now prepared to calculate the scalar perturbations in the synchronous- 
comoving gauge. We begin with the expression (6.19) for the scalar cf> 2 S . Again we split the 
quantities into a Newtonian and GR part, cj) 2 S = (/)+ </>(? s R . The only Newtonian part on the 
RHS of (6.19) arises from (3^ ps and some terms in T ij PS , where the latter is given in eq. (C.9). 
We thus have 
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(6.34) 


For the GR part we obtain after some tedious but straightforward calculations 
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To arrive at this result we have used the fact that 
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Ufom 0 
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n 2 n m0 


5P<7in + 3 


V 2 +E 


= Cu 


(6.36) 


is constant in time (and of course in space) which can be easily shown by taking the time 
derivative of the above and using the differential equations for the first-order and second-order 
growth functions, i.e., eqs. (5.12) and (5.31). A careful analysis reveals that Cm is vanishing. 

Now let us proceed with the derivation of E 2 S = E^ s + E^. From eq. (6.21) we obtain 
for the Newtonian part directly 
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To solve for the GR part, we obtain firstly 
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This expression can be drastically reduced by using the identity 

V- 2 V- 2 A,(^) = y + 2 ©0, 

and the relation (5.40). We then obtain 

^GR _ 20 

2s 9 u 2 n m0 

which concludes the derivation of E 2s . 


F>g m [(a n i - 1) + 3 @o] , 
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(6.38) 

(6.39) 

(6.40) 


Metric tensor: vector perturbations From the second-order gauge transformations (3.26) 
and (3.27) we obtain 

0 = uj 2ip + d 2ips , (6.41) 

F 2is = d 2lps + V- 2 T^ ps - v- 2 V- 2 T w Wipg , (6-42) 

where the solutions for ui 2 i p and T,; Jps are given in (5.42) and (C.9), respectively. For the 
first equation we have used the fact that Xj PS — V~ 2 £ fc fcips = 0, where £j PS is given in (C.6). 
The first equation gives then immediately the vector part of the gauge transformation 
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(6.43) 


where the transverse kernel TZi is given in (5.43). We have used (6.32) to solve for the time 
integral, and, as above, an integration constant was set to zero by requiring that at initial 
time the Eulerian and Lagrangian spatial coordinates coincide. Note that d 2 i PS is purely 
relativistic. 

Now we can easily derive from (6.42) the Newtonian and GR part of C 2 is=C 2 1 s + F 2 ™ 
They read 
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where we have used the definition ©o = —(1/2)V 2 V 2 D kl (ipofiipo,l), cf. with (5.21), and 
the relation (5.40). 


- 22 - 


























Metric tensor: tensor perturbations We finally derive the tensor modes in the metric 
from the transformation rule (3.28). As usual, we split X2ij s = X 2 \j s + X 2 ij s • an d obtain for 
the respective parts 
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(6.47) 


where 7r 2 jj P is given in Eq. (5.46), and we repeat the traceless and divergenceless tensor S 7J 
here for convenience 


SijiyX) — V To 5ij + T 2(<y?o,*j^ 74) PO,ik t P'o' t j') ■ (6.48) 

We shall comment on the tensors in section 8.2, when we summarise our results in the 
synchronous-comoving gauge. 

In concluding this section, let us emphasise that all the second-order synchronous- 
comoving gauge expressions obtained here are new. Only a few terms in these expressions 
were already known in the literature for the restricted case of an EdS universe (e.g., ref. [23]). 


Four-velocity As at first order, the scalar and transverse velocity vanishes in the La- 
grangian frame by definition. 

0 = v 2p - $ 2 PS + V- 2 O fc fcps , (6.49) 

0 = «4 P - 4s + ^ - V " 2 ^es • ( 6 ' 5 °) 

In both expressions, everything is already determined. We have used however these expres¬ 
sions to perform a consistency check of our calculations, e.g., from eq. (6.49) we can rederive 
v 2p (or, equivalently, by using the solutions of the RHS in (6.49) and (6.50), we can verify the 
validity of the comoving gauge conditions v = w l = 0 on the respective LHS’s, as it should 
for an irrotational and pressure-less fluid). For the perturbation of the zeroth component of 
the four velocity we obtain from eq. (2.7) that 

V 2 S = 0 • (6.51) 

Matter density Finally, from the formula (3.32), we find the second-order density contrast 
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(6.52) 

(6.53) 


The Newtonian part has been reported in the GR literature in e.g. [20], which however relies 
on a different solution technique namely the gradient expansion. The GR part, however, is 
fairly known see e.g., [6, 11, 13]. 
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7 From the synchronous-comoving gauge to the total matter gauge 

In this section we adopt the general formulas of section 3 to the transformation from the 
synchronous-comoving gauge to the total matter gauge [44] (sometimes also known as the 
velocity orthogonal gauge [27], or Eulerian gauge [6]). 

We define the total matter gauge perturbatively with (see section 4) 

E( r ) T = FT, = 0, (spatial gauge condition) (7-1) 

S l T° j( r ) T = 0 , (temporal gauge condition) (7-2) 

where S l = V~ 2 <9*. Relation (7.2) implies at first order (cf. eq. (4.6)) 

vi T + -Bi t = 0 , (73) 

and at second order (cf. eq. (4.8)) 

v 2t + B 2t - 2V~ 2 d l + 2(f>i T vi T j) = 0 . (7.4) 

As we shall show in the following, these conditions imply for an irrotational and pressure¬ 
less fluid that the time gauge generator of the gauge transformation from the synchronous- 
comoving gauge is vanishing at first and second order. 

7.1 First-order transformations 

Metric tensor From transformation rules derived in section 3.1, we have for the considered 
gauge transformation 


ipi T — T~Lcx i gT + di ST , (7f>) 

0i T = 0i s - ^V 2 /V , (7.6) 

Bi t = fii ST , (7.7) 

0 = Fi s + /?i ST . (7.8) 


Using the first-order result (6.11) for Fi s , equation (7.8) gives immediately for the spatial 
gauge generator 

a 2 7Vo 
PlST “ 3 n 2 n m0 ■ 

Note that, apart from the minus sign, this spatial generator is the same as from the Poisson 
gauge to the synchronous-comoving. The minus sign arises because in this section we perform 
the transformation from Lagrangian to Eulerian frame, whereas in section 5 it was the other 
way around. 

To obtain the temporal gauge generator ai ST , we sum up eqs. (3.13) and (3.19) which 
gives the general expression vi + Bi = v\ + B\ — aq. For the considered gauge transformation, 
where iq g = B i g = 0 and because of the temporal gauge condition (7.3), this implies 



cti gT — 0. 


(7.10) 
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Since we assume vanishing first-order vector perturbations, we have d\ ST = 0. Having deter¬ 
mined the gauge generators, the remaining first-order scalar perturbations are easily derived, 


Vh T = 0, 

/ _ 5 

tIt — 2 *Pin ) 

_ 2 Vtp 0 
lT ~ 3 n%n m0 ' 


(7.11) 

(7.12) 

(7.13) 


Four-velocity 

by 


The first-order spatial peculiar velocity in the total matter gauge is given 


^1t /^lST 


2 7>(/?q 

3ft§n m0 ’ 


(7.14) 


which follows directly from (3.19). Note that the spatial peculiar velocity (7.14) is identical 
to the one in the Poisson gauge (cf. (5.16)). For the perturbation of the zeroth component of 
the four velocity, since ipi T = 0 we obtain from eq. (2.7) that 


v 


0 

1 T 


= 0. 


(7.15) 


Matter density The matter density is identical to the one in the synchronous-comoving 
gauge because the time coordinate is identical and the matter density is unaffected by a purely 
spatial transformation at first order. We thus have directly from (3.21): 


<5i t 


_ 2 PV 2 yp 
ls “ 3?7pW 


(7.16) 


7.2 Second-order transformations 

Metric tensor: scalar perturbations From the transformation rules at second order 


(3.22)-(3.25) we have 

ip2 T = 'Ha2 ST + d2 ST + n sT , (7.17) 

= ^2 S - g T \. ST --V 2 /3 2st , (7.18) 

H 2T =/3 2sT + V- 2 E fc ifcsT , (7.19) 

0 = E 2s + /3 2st + ^ V- 2 V- 2 Tf T .. - ^ V~ 2 T\. st , (7.20) 


where the first-order squared terms Hst, T?j st , and Sj gT are given in appendix C. In writing 
down these equations we have only used the spatial gauge condition for the total matter 
gauge, see eq. (7.1). The temporal gauge condition (7.4) of the total matter gauge becomes 
at second order 


v 2t + B 2t = 2V 2 d l (^iT-Bix.t + 2<f>! T v lTi i) = . (7-21) 

y /T-g'^mO 

where we have used the first-order results from section 7.1. Now, to obtain the temporal 
gauge generator a 2sT , we sum up eqs. (3.24) and (3.30) which yields the general expression 
v 2 + B 2 = v 2 + B 2 + V _2 H fc . fc — a 2 + V 2 S fc ; fc. For the considered gauge transformation, 
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where V 2 S = i? 2 s = f2g T = 0 and V 2 Eg T fe = — 9 ^ 2 ^ in o y?o (see appendix C), we obtain for 
the time gauge generator at second order 


Ct2oT — ^ ■ 


(7.22) 


Also the derivation of the longitudinal part of the spatial gauge generator is straightforward. 
From eq. (7.20) we get 
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(7.23) 

(7.24) 


where /? 2 ST = ■ The Newtonian part is, of course apart from a global minus sign, 

identical with the Newtonian part of /32 PS , see eq. (6.30). The GR part is identical to 
(see eq. (6.31)), apart from the terms ~ T) 2 , which arise from the different time coordinate 
used. 

The calculation of the other metric perturbations is now trivial. We obtain 
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(7.29) 


9 n 2 n m0 /ru 377 2 54 m0 

Note that eq. (7.28) is the second-order peculiar velocity potential in NEPT (see eq. (5.50)). 

Metric tensor: vector perturbations From the second-order gauge transformations (3.26) 
and (3.27) we obtain 
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(7.30) 

(7.31) 


where Tjj ST and Ej ST are given in eqs. (C.12) and (C.14). Again, the calculations are simple 
so we leave them as an exercise. We obtain 
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where is given in eq. (5.43). 
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Metric tensor: tensor perturbations We obtain this simple result 

X2ijT = + ’ (7 ' 34) 

where S t] is given in eq. (6.48), and ir 2 ij P , which contain propagating tensor modes (gravita¬ 
tional waves) in the Poisson gauge, can be found in (5.46). Observe that the “Newtonian-like” 
tensors (growing with the same amplitude as second-order Newtonian perturbations), and 
partly the post-Newtonian tensor, that we have in the synchronous-comoving gauge, disap¬ 
peared. 

Let us emphasise again that all these ACDM results are new. 


Four-velocity The transformation rules (3.30)-(3.31) deliver simple expressions, for the 
scalar and vector part of the spatial peculiar velocity we find firstly 

V2 T = “/Wr > (7-35) 
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(7.37) 

(7.38) 

(7.39) 


where we split as usual V 2 T = + v^. For the perturbation of the zeroth component of 

the four velocity we obtain from eq. (2.7) that 




= 0. 


(7.40) 


Matter density Finally, according to (3.32) the matter density perturbation transforms as 


- h s + 2<5i S) fc£i gT • 


(7.41) 


Since the last term on the RHS is purely Newtonian, the GR corrections are identical to the 
one in the synchronous-comoving gauge. 11 We get then for the density contrast 
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(7.42) 


which agrees with the findings in [6]. 

14 We note that it is in principle possible to define an Eulerian gauge where the last term on the RHS 
of (7.41) will not only contain a Newtonian part but also a relativistic part [47]. 
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8 Summary of the results 


In this section we summarise our results in the three gauges, namely the Poisson gauge 8.1, 
the synchronous-comoving gauge 8.2 and the total matter gauge 8.3. Here we only show the 
results for the metric and for the density contrast. For the summary of the gauge generators 
see section 8.4, and for more results we refer to the sections above. For an overview of the 
used notation, see tab. 1 in appendix A. We recall that for the EdS limit (H m = 1) we have 
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V = P , 
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(EdS universe) (8.1) 


8.1 Poisson gauge 

We obtain for the metric components in the Poisson gauge 
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where TT 2 ij P is given in eq. (5.46). The density contrast is 
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8.2 Synchronous-comoving gauge 


We obtain for the metric components in the synchronous-comoving gauge 
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where the tensor 
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(8.9) 


is the solution of the gravitational wave equation 


7T 2ijs + ‘IT-Lit 2ij s 
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( 8 . 10 ) 


This result, valid for ACDM, is new and coincides with ref. [23] in the EdS limit. Comparing 
the gravitational wave equation (8.10) with the one obtained in the Poisson gauge (Eq. (5.44)), 
it is evident that the nature of the wave equation has not changed but only its source term. 
The solution of the wave equation includes Tt 2 ij P , given in eq. (5.46), however, as a consequence 
of the different source term, solution (8.9) contains also additional tensor perturbations that 
grow with J 7 , V and V 2 . The physical interpretation of the various tensor perturbations is 
highly non-trivial but has been attempted in refs. [6, 23]. 

For the density contrast we obtain 
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8.3 Total matter gauge 

We obtain for the metric components in the total matter gauge 
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For the density contrast we obtain 
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8.4 The gauge generators 

From the Poisson gauge to the synchronous-comoving gauge We use (3.2) and then 
obtain the fluid trajectory as viewed from an observer in the Poisson gauge 
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where (r, q) are the Lagrangian coordinates. The relation between the temporal coordinates 
is 
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From the synchronous-comoving gauge to the total matter gauge We use (3.1) and 
then obtain the fluid trajectory as viewed from an observer in the total matter gauge 
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For the temporal coordinate we simply have = t. 


(8.18) 
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Comparison with the Newtonian trajectory In Newtonian gravity, the Newtonian 
trajectory is up to second order 
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where we recall that bars denote Euclidean coordinates. 


(8.19) 


9 Conclusions 

We have considered relativistic perturbations for a collision-less and irrotational fluid up to 
second order in a ACDM Universe. In detail, we have determined the fastest growing mode 
solutions of all metric perturbations, the density and the velocity for three common gauges, 
namely for the Poisson gauge, the total matter gauge and for the synchronous-comoving gauge. 
We considered growing mode initial conditions with primordial non-Gaussianity. First, we 
have confirmed the metric expressions of [41] in the Poisson gauge and corrected for some 
typos. We found major simplifications for the expressions in the Poisson gauge. Let us 
emphasise that, because of these findings, all perturbations up to second order are fully 
described by only two time coefficients, which are identical with the first-order and the second- 
order growth of the Newtonian displacement held. Furthermore, our resulting solutions are 
very compact and come with simple physical interpretations. 

Having determined the solutions in the Poisson gauge, we have performed a gauge trans¬ 
formation to the synchronous-comoving gauge. This gauge transformation can be understood, 
in terms of fluid mechanics, as a perturbative coordinate transformation from a specific Eule- 
rian frame to the unique Lagrangian frame. Our results in the synchronous-comoving gauge 
are, to our knowledge, entirely new, only the second-order density was known so far in the 
literature [11, 24-26]. 

Then, we have considered the gauge transformation from the synchronous-comoving 
gauge to the total matter gauge, which, as above, can be understood in terms of a coordinate 
transformation, here from the unique Lagrangian to another Eulerian frame. This particular 
Eulerian frame is very important, as the resulting time coordinate in the total matter gauge 
is identical with the one in the synchronous-comoving gauge. That is, this particular Eulerian 
frame makes use of the proper time of the fluid particles. This choice of Eulerian frame is thus 
very convenient when relating GR results to Newtonian investigations, since the relativistic 
Eulerian-Lagrangian correspondence makes use of the same time coordinate. So to say, the 
problem of relating a GR description to a Newtonian world has been reduced from a 4D 
problem to a 3D problem. Such considerations are for example useful when generating GR 
initial conditions for Newtonian IV-body simulations (see [47]), the latter being one area of 
application of our paper. 

Our results for the coordinate transformations in section 8.4 dictate the relativistic tra¬ 
jectories of CDM particles in a ACDM Universe, and these results tell us explicitly how the 
Newtonian particle trajectories are deformed because of general relativistic effects. These GR 
corrections could be incorporated into N-body solvers with the aim to include GR effects in 
Newtonian N-body simulations effectively up to second order in relativistic perturbation the¬ 
ory. Our results for the density and velocity in three common gauge choices have also direct 
outcomes for relating theory with observations of the large scale structure of the Universe. 
Specifically, measurements of the galaxy number density require both the knowledge of gauge 
effects (stemming from the choice of time-like hypersurfaces) and of GR effects due the fact 
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that observations are performed on the past light-cone (e.g., observed photons are gravita¬ 
tionally distorted while passing through a clumpy Universe), the former being performed in 
the current paper and the latter being investigated elsewhere (see e.g., [48, 49]). 

Finally, we have no reason to hide the pleasure that all our results are fairly simple. Of 
course, the reason for such simple expressions is hidden in the laborious work of solving the 
time differential equations in the Poisson gauge, which were not known so far. 
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A Used notation 

In tab. 1 we give an overview of the notation used in this work. For the case of an EdS 
universe, we have f2 m = 1, a = rj 2 , %/2r/, %o^mO = 4, g = g m = 1, V = rj 2 , and T = (3/7 )r/ 4 . 

B Lie derivatives 

The Lie derivatives and £| along the vector held have the following expressions: 


• scalars 


A/ = f,uC 

c\f = + C + fA a 


(B.l) 

(B.2) 


• vectors 



(B.3) 

(B.4) 


• 1-forms 



(B.5) 

(B.6) 


• metric tensor 



(B.7) 


(B.8) 
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metric (/j, = 0,1,2, 3) 

eq. (2.1) 

9ij 

spatial metric (i = 1, 2,3) 

eq. (2.5) 

a (v) 

cosmic scale factor 

eq. (2.2) 

Ti = a/a 

conformal Hubble parameter 

eq. (5.3) 

Om0?) 

matter density 

eq. (5.6) 

V( V ) 

structure growth of linear density fluctuations 

eq. (5.13) 

Hv) 

time coefficient of second-order Newtonian displacement 

eq. (5.31) 

g(v) 

growth suppression rate (= P/a) 

eq. (5.11) 

f = V/(HV) 

differential growth 

eq. (5.19) 

S(v,x) 

density contrast (5 = (p — p)/p) 

eq. (2.9) 

= g(v)<Po 

linear gravitational potential at time r/ 

eq. (5.9) 

<Po(vo,x) 

tp at the present time rjo 

eq. (5.9) 

74n = gin^O 

linear gravitational potential at initial time 

eq. (5.1) 

an 1 = (3/5)/„! + 1 

primordial local non-Gaussianity 

eq. (5.1) 

v 2 

inverse spatial Laplacian with Euclidean metric 

eq. (3.28) 

Dij 

longitudinal extraction operator (Dij = didj — (1/3)V 2 ^) 

eq. (2.4) 

^o(jlQ,x) 

second-order “Newtonian” kernel for the displacement 

eq. (5.38) 

@o(vo,x) 

second-order “GR” kernel (V 2 ©o = 1 J / o — (l/3)(V</?o) 2 )) 

eq. (5.21) 

TZi(r] 0 ,x ) 

second-order “GR” vector kernel 

eq. (5.43) 

Sij{r) o,x) 

second-order “GR” tensor kernel 

eq. (5.45) 

K2ij P (tj,x) 

secondary tensor perturbations in the Poisson gauge 

eq. (5.46) 

i\> 

scalar perturbation in <700 

eq. (2.1) 

B 

scalar perturbation in < 70 * 

eq. (2.3) 

Ui 

vector perturbation in goi 

eq. (2.3) 

<t> 

scalar perturbation in the trace of < 7 ^ 

eq. (2.4) 

E 

scalar perturbation in the trace-less part of gij 

eq. (2.4) 

Ft 

vector perturbation in gij 

eq. (2.4) 

Xij 

tensor perturbation in g^ 

eq. (2.4) 

ufl = K + v^)/a 

4-velocity (u fJ 'u u g fll/ = — 1) with perturbation v ^ 

eq. (2.6) 

v o 

perturbation in the time-component of u^ 

eq. (2.7) 

Vi = V }i + Wi 

(longitudinal and transverse) spatial velocity perturbation 

eq. (2.8) 

? = (?,?) 

gauge generator 

eq. (3.1) 

P 

II 

i/n 

o 

time component of the gauge generator 

eq. (3.3) 

+ 

II 

Uy 

(longitudinal and transverse part of) spatial gauge generator 

eq. (3.3) 


Table 1. Used notation in this work. 
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C First-order squared terms arising at second order 


The contributions from products of first-order terms to the second-order gauge transforma¬ 
tions in section 3.2 are given by 


n = ai 


a i + hUai + (H + 2'H 2 )a\ + AT-Lipi + 2-01 


+ 2di(di + 2'0i) + Ci (oq + Hoti + 2ip\) k + Ci oi\,k ~ 2,Bi k — Cifc 


S, = 2 


( 21-LBii + Bu)ai + Bu^kii ~ 2Vh«i,i + -BifcCgi + Bu&i + 2Ciik£i 


+ AUaidu - a M ) + di(6i - Sa^i) + ai(Cn - «i ,i) 

+ Ci (Ci i,k + 2Ci k,i) + Ci (Cii.fc “ - ai.fcCi,* > 

= Cu (2^1 + «i + 2Tiai) — aiCii — CiCii.fc + CiCii.fc 

- 2ai (v u + Uvu) + 2w lijfc ^ - 2v^ li>k , 


(C.l) 


(C.2) 

(C.3) 


and 


T ij — 2 


(H + 2'H 2 )a\ + %(aidi + ai^Ci) 


Jij 


— 2a\^a\j + 4ai(Cijj + 2T-LCuj) + ABi^a^ 


+ 4 


CiijjfcCi + 2C'ifc(iCij) + 8"HaiCi(ij) + 2Cife,jCij 


+ 2aiCi(jj) + 2Ci(ij)fcCi + 2Ci(i,fcCij) + 2Ci(iCtij) 


(C.4) 


In the next two subsections we give the explicit expressions for our specific transformations. 


From the Poisson to the synchronous-comoving gauge The second-order contribu¬ 
tions which arrive solely from first-order squared terms are in the case Poisson gauge to 












synchronous-comoving gauge 


n PS = aipsi'^o - 2^Vo + /^lpgfi'^o ,k = ~ aip S ^o + \vin9Vo ~ 3c/Vo + Pi* s 9Vo,k 

CL O 

2P 2 9 \ , 2P 2 


I £ Ls 2 1 2 £ Ls l 

— I q 99in — Q nj2n ) V^O Q n/2n ^0»^b 

\ 3 3aT~L^il m Q j 3tt7^QJl m o 


(C.5) 


k 

£*PS 6^(/9o^ips,^ “I” 27/ai PS ai PSj i + 2o'i PS} / Ci 5^p S ^ 

^S^0^1ps,£ g V^in^lps,^ “I - ^^lpsj^^lps,i 

8 £>X> 20 £ in X> 8 VV k 

= 3o^fl m0 WW + i^ WOii + 9 ('H 2 n m0 ) 2ipo ’ kip ’°’ i 

5 

H^PS ^^O^lpg,! 4~ 2^QJipgC^lpQji “1“ 3^^o,i 

_ 20 £ in P . PP 


(C.6) 


;®lpS — g V^inC^lpg,* 5cCl p g “I - 6ttl 


-PS 74* 


20 C/i n P . VV 

~9 ~ 4 anfamo* 0 * 0 ’* ’ 


(C.7) 


T — 

- 1 *?PS 


5ij 


2^a? PS - 6Uai ps g<p 0 + 2%ai ps , fc /3ip g - 4ai PS ^(/? 0 - 4gv?o,fc/3£ s 

A: 

— 8<7<y9o/?i P g4~ 87-toipg/3i P g,ij 4“ ^/^ips,fc*/^i P gj 4” 2oipgCn PS ,jj 4~ 2/3\ps,ijkfi\ FS 


20 g m V )fc ,4 P 2 ,* 

9 'H 2 n m0 ^ kip ° 3a^gn m0 V0 ^° 


100 O 10 „ 2 4 P 2 \ 2 20 C/inP 

- jm« - 29 + 3^2^] *>„ + 

80 S1 „P 8 P 2 16 P 2 t 8 P 2 

+ 9(WpUp W ”-‘ J + ¥(WM Vo ' il ^ + 


(C.8) 

k 


°ij 


9 (H 2 n m0 ) 2 ^°’ feii ^ 
(C.9) 

To derive the second equality of each expression we have made use of eq. (6.13). It is worth¬ 
while to derive the trace of Ty PS explicitly: 

7 100 9 9 X 

-g-0in - 10 5.9iti - 6 g 2 + 4- 


'Y'k 

1 fc PS — 


aH 2 Q m0 J ^ 


20 ffmP 


-w k P^ A: 30 3inP _2 

+ TT2?r^^ + 4 -i,2 0 yo.fcyb + -qyJTj ^° V 

o rtgJimO ®'t|-|Si m o J rtnSimO 


tg“m0 
P 2 


tg^mO 
16 P 2 


tg' hnO 

„h + 8. ® 2 


8 P 2 16 P 2 u 8 P 2 t , 

+ 9(ii|0rf f,T ^ + ^ ' (CA0) 

From the synchronous-comoving to the total matter gauge Since there is no time 
transformation involved, the second-order contributions from first-order squared terms reduce 
to 


TT =—P k P =—— ^ 

1A ST ?lsT^ lfc ST 0 " 


"9(H 2 52 m0 )2^ 0 ’ fc ^ 

„ \ Pk Pk (/- r\/- \ /-k P 40 'VQ\ j3 

^*ST If^lifcsslgx 4“ Slg T (4l«ST,fc 4” 2£lfcg T ,i) 4" Slg T ?lig T ,fc ^ q ^Po^Po,i 

H*ST 9 , 


(C.ll) 

(C.12) 

(C.13) 
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and 


T =4 

*?ST ^ 


C'lyS.fcClsx + 2C lk(i S €l ST ,j) + 2 £lfcST,i£l S T,j + 2 ^lsT(*,i)fc^l + 2 ^lsT(j,fe^l S Tj) 


40 T)g m ^ 80 V 9] 

6ij ~ 


8 V 


jvj is iym o x--' / _L , 

~ 9 (?^Q m o) 2 ( v 2^ 0 ’ fcV9 ° + 


„,fe 


(C.14) 


where for the second equality we have substituted the various expressions from the metric in 
the synchronous-comoving gauge and from the first-order transformations. It is worthwile to 
derive the trace of Yjj ST explicitly: 


'Y'fc 
1 ksT 


40 


Yh/in k 

PO,kPo 


3 n 2 n m0 


80 Y>g in 

9 HZSlmD 


^oVVo - x 


V 2 


9 (uln 


mO] 


-^ 2 (<Po,k<p{) + l Po,ki l Po l 


(C.15) 
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